Heavy-fermion spin liquid in the strong hybridization limit of the finite-U Anderson 
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Studying the finite-U Anderson lattice model in the strong hybridization limit, we find a heavy- 
fermion spin liquid phase, where both conduction and localized fermions are strongly hybridized 
to form heavy fermions but this heavy-fermion phase corresponds to a symmetric Mott insulating 
state owing to the presence of charge gap, resulting from large Hubbard-U interactions in localized 
fermions. We show that this heavy-fermion spin liquid phase differs from the "fractionalized" Fermi 
liquid state, where the latter corresponds to a metallic state with a small Fermi surface of conduction 
electrons while localized fermions decouple from conduction electrons to form a spin liquid state. 
We discuss the stability of this anomalous spin liquid phase against antiferromagnetic ordering and 
gauge fluctuations, in particular, instanton effects associated with confinement of slave particles. 
Furthermore, we propose a variational wave function to check its existence from the microscopic 
model. 
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I. INTRODUCTION 

Spin liquid has been intensively studied, motivated 
from the resonating-valance-bond (RVB) scenario for 
anomalous finite temperature physics of the Pseudogap 
phase in high T c cuprates.[l] However, such spin liq- 
uid physics is fascinating itself because this phase has 
a nontrivial order called topological or quantum order 
beyond the description of the Landau-Ginzburg- Wilson 
paradigm for phase transitions, and its low energy physics 
is described by a gauge theory, far from the conventional 
structure of condensed matter theories, allowing spin- 
fractionalized excitations called spinons.[2] 

Recently, several geometrically frustrated insulators 
are proposed to be genuine symmetric Mott insulating 
phases, i.e., spin liquids. [3] Although such magnetic frus- 
tration can be a strong candidate as the mechanism for 
the existence of spin liquid, it is still interesting to find 
another mechanism for its existence. In this paper we 
propose a heavy-fermion spin liquid phase in the strong 
hybridization limit of the finite-U Anderson lattice model 
(ALM). 

An important thing is how to suppress magnetic or- 
dering. Antiferromagnetic order is well known to occur 
at half filling in the square lattice due to Fermi-nesting. 
One way to suppress the antiferromagnetic ordering is to 
make spin singlets introducing another band-electrons. 
Furthermore, if the nesting property can be destroyed, 
such an ordering tendency will be much more suppressed. 
This motivates us to consider the ALM since it shows a 
large Fermi surface when conduction electrons are hy- 
bridized with localized electrons. 

In this paper we consider the strong hybridization 
limit, turning on on-site Hubbard interactions to local- 
ized electrons. Two possibilities are expected. One is 
that strong local interactions weaken hybridization, caus- 
ing the "fractionalized" Fermi liquid phase, where local- 
ized fermions form a spin liquid phase while conduction 



electrons are in the Fermi liquid state. [4] However, Fermi- 
nesting in the localized fermion-band would result in an- 
tiferromagnetic ordering. Magnetically frustrated inter- 
actions are required for the fractionalized Fermi liquid 
to realize. The other is that the hybridization still sur- 
vives against local interactions, but charge fluctuations 
are suppressed due to such interactions. Since Fermi- 
nesting does not exist and Kondo singlets are formed 
owing to hybridization, antiferromagnetic ordering does 
not arise, but this corresponds to an insulating state ow- 
ing to charge gap. Because such a phase contains all 
symmetries of the original lattice model, it can be identi- 
fied with a spin liquid state with heavy neutral fermions, 
i.e., spinons. We propose that this heavy-fermion spin 
liquid state can emerge in the strong hybridization limit 
of the finite-U ALM. 

To describe Mott transition, it is necessary to take 
on-site Hubbard interactions non-perturbatively. Re- 
cent studies on the Hubbard model have revealed that 
the slave-rotor representation can treat such interac- 
tions well, describing the Mott transition. [5] Applying 
the U(l) slave-rotor representation to the finite-U ALM, 
and performing gauge transformation [6] appropriate in 
the strong hybridization limit, we find an effective U(l) 
gauge Lagrangian in terms of renormalized conduction 
and localized fermions and collective density-fluctuation 
bosons interacting via U(l) slave-rotor gauge fluctua- 
tions. We show that such density fluctuations associated 
with localized fermions become gapped, increasing the 
Hubbard interaction-U in the strong hybridization limit. 
This gapped phase is identified with an insulating phase 
with all symmetries of the ALM, thus heavy-fermion spin 
liquid. We discuss the stability of the spin liquid phase 
against gauge fluctuations, in particular, instanton ex- 
citations associated with confinement of slave-particles. 
Furthermore, we propose a variational wave function to 
check its existence from not the effective field theory but 
the microscopic model itself. 
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It should be noted that the heavy-fermion spin liq- 
uid phase differs from the fractionalizcd Fermi liquid 
state [4] because the former is an insulator, thus there 
is no Fermi surface while the latter is a metal to have 
a small Fermi surface of conduction electrons. In this 
respect the present Mott transition in the ALM should 
be discriminated from Kondo breakdown as the orbital- 
selective Mott transition[7]. Furthermore, the spin liquid 
state in the fractionalizcd Fermi liquid phase would be 
unstable against antiferromagnetic ordering owing to the 
presence of Fermi-nesting. 

It is important to notice that the limit considered in 
this paper is different from that studied in the context of 
heavy fermion physics. [8, 9] Quantum phase transitions 
in heavy fermion compounds are studied in the large-U 
limit of the ALM, where the Kondo lattice model (KLM) 
or infinite-U ALM is in the main interest. On the other 
hand, the present study considers the strong hybridiza- 
tion limit, not seriously taken into account in the KLM 
or infinitc-U ALM context. 



II. U(l) SLAVE-ROTOR REPRESENTATION OF 
THE FINITE-U ANDERSON LATTICE MODEL 

We start from the finite-U ALM 

LALM = J2 C L( d r ~ - t E (4a C 3* + H - C ) 

-Vj2(cld t « + H.c.) + J2 4a(dr + e d )d ia 

ia ia 

+Uj2 d hd l1 4 l d ll , (1) 

i 

where Cj CT represents conduction electrons with dispersion 
e c k = —2t(cosk x + cos k y ) in two dimensions, and di a 
localized electrons with localized level e d . 

Quantum phase transitions are expected to occur via 
competition between the hybridization term with V and 
the local interaction term with U. Without the interac- 
tion term of U this model is exactly solvable, resulting 
in hybridization between the conduction and localized 
bands. Our problem is to study what happens to the 
hybridization when the local interaction U is turned on. 

The local interaction term can be decomposed into the 
charge and spin channels 

uj24^4idn = |EE! d U„ - !) 2 

i i a 

-t BE ad ^ 2 + ? E d - d - - t E ^ 

i a icr i 

In this paper we do not consider the spin channel, jus- 
tified in the strong hybridization limit where antiferro- 
magnetic ordering is severely suppressed. 

Recently, we could describe the genuine Mott transi- 
tion without symmetry breaking, introducing the slave- 
rotor representation in order to take into account local 



interactions non-perturbatively in the Hubbard model. [5] 
Decomposing the localized electron di a — e. i f]i a , one 
can rewrite the ALM in the following way 

Lalm = E 4x(<9t - M)c M - t E (4<7 C jt + H - c -) 

ia (ij)<j 

-V EKt^N- + H.c.) + E vl(dr + e d ) VilT 



+tE l ?-*E l ^ 

i i 

+i^ipi(Li - vtviv - 1]), 



(3) 



where e d is replaced with e d — U/2. It is easy to show 
that Eq. (3) is exactly the same as Eq. (1) with Eq. (2) 
after integrating out the ipi field with rj ia = e l0i d ia . In 
this expression the electron Hilbert space |d i(T > is given 
by the direct product of the fermion and boson Hilbert 
spaces \r)i a > ® |Lj > according to the decomposition 
die = e~ %6i r]ia, where Li represents an electron density at 
site i. It is clear that any decomposition method enlarges 
the original electron Hilbert space, thus an appropriate 
constraint associated with the decomposition should be 
imposed. The Lagrange multiplier field tpi expresses the 
U(l) slave-rotor constraint Li = J2 a ViaViv — 1, implying 
that the fermion and boson Hilbert spaces are not inde- 
pendent, thus the two operators r\i G and e lSi also. Then, 
e -«£>i j s identified with an annihilation operator of an elec- 
tron charge owing to the constraint Li = J2 a vj^Via — 1 
and the canonical relation [Li, Oj] = —iSij imposed by 
—iLid T 6i. In this respect collective density fluctuations 
associated with localized fermions can be taken into ac- 
count non-perturbatively in the U(l) slave-rotor repre- 
sentation. 



III. HEAVY FERMION PHYSICS: KONDO 
BREAKDOWN AS THE ORBITAL-SELECTIVE 
MOTT TRANSITION 

The problem is how to treat the hybridization term 



with 



,-iOt 



One direct way is to integrate out con- 



duction electrons and obtain an effective Lagrangian in 
terms of localized fermions and their density-fluctuation 
bosons. [4] Integrating out the density field Li and con- 
duction electron field Ci a , we obtain 

L eff = E ^iraidr + - » E V* E ~ ^ 

ia i a 

i 

(4) 



i j a 



where A? ■ TT , is the single particle propagator of the con- 
duction electron, given by 

V 2 



iuj + n - e c 
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in the energy-momentum space. 

Although the above treatment itself is exact so far, it 
has an important assumption that the identity of conduc- 
tion electrons is sustained. In other words, feedback ef- 
fects of localized fermions and density fluctuations to con- 
duction electrons, self-energy corrections of conduction 
electrons, are not introduced. Since such feedback effects 
can be induced only via the hybridization coupling term, 
this treatment may be regarded as the weak hybridiza- 
tion approach, not justified in the strong hybridization 
limit. [10] Alternatively, it may be viewed as the large- U- 
limit approach owing to small V/U, where the infinite- U 
ALM can be considered or the KLM is taken via virtual 
charge fluctuations. This type of treatment has been 
utilized both intensively and extensively, associated with 
the single-impurity [11] and heavy fermion physics[9]. 

Performing the Hubbard-Stratonovich transformation 
for the non-local hopping term in both time and space, 
we obtain an effective Lagrangian 

L ?ff = E ipr + td)m* - « E ^ E vhvi* - 1] 

icr i cr 

+ E < (°r - m + e c q ) x 6 q + \ E^ - Vl f 

q i 

-I'EEix'v + e'^e-^], (5) 

ij cr 

where x\ an d X% are effective hopping parameters de- 
termined in the self-consistent analysis 



X% = V(e»'{d T -» + <e j )- 1 e-»'), 

4 = ^t(^-/j+4rv>- 



(6) 



It is not easy to perform the self-consistent analysis 
because such effective hopping parameters have both fre- 
quency and momentum dependencies (nonlocal in time 
and space). This non- locality may give rise to crucial ef- 
fects to the heavy fermion physics. Since this is not our 
main subject, we leave its detailed analysis in an impor- 
tant future problem, and here, discuss what is expected 
in this treatment. 

Generally speaking, two phases would be allowed in 
this approximation, determined by the condensation 
of slave- rotor bosons. In the large V/U limit collec- 
tive density-fluctuation bosons become condensed, and 
heavy-fermion Fermi liquid appears to form a large Fermi 
surface since (e l9i ) ^ results in (cj^Tftcr) ^ as shown 
in the hybridization term of Eq. (3). In the small V/U 
limit such boson excitations become gapped, and local- 
ized fermions decouple from conduction electrons owing 
to (c\ a r)i a ) = 0, forming a spin liquid state. Such a phase 
is called the fractionalized Fermi liquid state with a small 
Fermi surface of conduction electrons. [4] As a result, the 
heavy-fermion Fermi liquid to fractionalized Fermi liq- 
uid transition is identified with Kondo breakdown as the 
orbital-selective Mott transition [7] in the U(l) slave-rotor 
representation of the fmite-U ALM. However, the spin 



liquid state would be unstable against antiferromagnetic 
ordering owing to Fcrmi-ncsting in the square lattice. 

Although one important issue, that is, the abrupt 
volume-change of the Fermi surface at the heavy-fermion 
quantum critical point, [12] can be explained by the 
Kondo breakdown transition as discussed above, the 
fact that the antiferromagnetic transition of localized 
fermions arises simultaneously at the same quantum crit- 
ical point as the Kondo breakdown transition is still far 
from our understanding because such a phenomenon is 
beyond the description of the Landau-Ginzburg- Wilson 
theoretical framework for phase transitions. [13] How to 
incorporate such an antiferromagnetic transition in the 
Kondo breakdown one is an important open problem. 



IV. STRONG HYBRIDIZATION APPROACH 

In the previous section we have discussed the heavy 
fermion physics in the "large-U" (more accurately, weak 
hybridization) limit of the ALM although its mathemat- 
ical construction is performed at fmite-U. In this section 
we take the strong hybridization approach, considering 
that the hybridization term is relevant in the renormal- 
ization group sense when U = 0. We show that this 
type of approach allows a new possible phase called the 
heavy-fermion spin liquid in the finite-U ALM. 



A. Gauge transformation and effective Lagrangian 

The relevance of the hybridization coupling in U = 
motivates us to consider the gauge transformation 
of Ci a = e~ l9i ^ ia , where this type of approach has 
been utilized in the context of quantum disordered 
superconductors[14]. Integrating out the density field Li 
and inserting another slave-rotor decomposition Ci a = 
e~ %di ^ia into Eq. (3), we find the following expression 



Dip icr Dr] icr D9iD(pi exp 



-H - id r 9i)^ ia - t ]T (V-Le^e-^Vv + H.c.) 



icr icr 

-iE^E^U- - ^ + jjY.^ w) 2 }]- (7) 



Couplings between charge fluctuations and renormal- 
ized conduction electrons can be decomposed using the 
Hubbard-Stratonovich transformation. Following Ref. 



4 



[5] , we find the effective Lagrangian from Eq. (7) 
L e ff = Lq + L^p + L n + Lv + L$, 

Lo = j-J E Q ri ~ 7? E Qri(tPri ~ Pi) + ^ 



U ^ ™ U 

i 

£ V> = E ^t( d r - V)4>i* ~ i E Pi E - !] 

ia i <j 

ia i a 

L v = -Vj^WUi* + H.c), 

ia 

L 9 = jj ^2(d T 0i - ¥n) 2 - 2ty^2 / cos{9 i - 6j + ay), 



where the hopping parameters are represented as Xij = 
xe laij and yy = ye WiJ , and ipi = <p ri — q ri is used with 
q r i = Uqi/2. The unidentified parameters can be deter- 
mined self-consistcntly in the saddle-point analysis 

Pi = (d T 6i), q ri = i-^(^2i>LA<r - 1), 

a 

x= \(e- ie <e»i+H.c.)\, y=\(^Uia + H.c.)\. 

a 

(9) 

Pi and q r i play the role of renormalization for the fermion 
and boson chemical potentials, respectively. 

The effective Lagrangian Eq. (8) is the main result of 
this paper, showing an important feature of the finite-U 
ALM, that is, the presence of two kinds of relevant inter- 
actions. As mentioned before, the hybridization coupling 
is relevant to give rise to a band-hybridized metal if the 
local charge-fluctuation energy is not taken into account. 
On the other hand, an appropriate treatment of local 
charge fluctuations has been shown to result in the Mott 
transition from a spin liquid Mott insulator to a Fermi 
liquid metal at a critical value U c without hybridization, 
i.e., in the Hubbard model. [5] 

Compared to the effective Lagrangian Eq. (5) for the 
heavy-fermion quantum phase transition, the effective 
Lagrangian Eq. (8) always allows the band-hybridization 
while Eq. (5) does not. In Eq. (5) the boson con- 
densation corresponds to the hybridization transition 
since it gives an effective hybridization coupling constant 
Veff — \(e~ l9i )\V. Such a transition is involved with 
the abrupt volume change of the Fermi surface. On the 
other hand, in Eq. (8) the renormalized conduction and 
localized fermions are always hybridized, forming a large 
Fermi surface of spinons (neutral fermions). Then, the 
boson condensation in Eq. (8) corresponds to the genuine 
Mott transition in these heavy fermions. 

Since the slave-rotor variable has its own dynamics in 
the saddle-point approximation, given by the U(l) ro- 
tor (XY) model, the quantum transition described by 



condensation of the rotor field occurs when Dy/U ~ 1 
with the half bandwidth D, where Dy is an effective 
bandwidth for the rotor variable. [15] It is important to 
understand that the ratio V/U between the hybridiza- 
tion and on-site interactions controls the hopping pa- 
rameter y, where the hopping parameter decreases as 
this ratio decreases. This results in the quantum crit- 
ical point (V/U) c between the heavy-fermion Fermi liq- 
uid and heavy-fermion spin liquid, which differs from the 
heavy-fermion quantum critical point between the heavy- 
fermion and fractionalized Fermi liquids. 

^' It is interesting to notice that this treatment [Eq. (8)] 
gives rise to two different energy scales. [16] One corre- 
sponds to the band- hybridization scale T K , and the other 
is associated with the coherence scale T FL , resulting in 

.o\Fermi liquid. Tk is proportional to the hybridization gap 
V in the strong coupling approach, consistent with the 
single impurity energy scale in the strong coupling limit. 
On the other hand, Tfl is proportional to the effective 
stiffness parameter Dy\ (e l0i ) | 2 controlling the phase co- 
herence of [17] 



B. Mean-field analysis 

For the saddle-point analysis we resort to large N gen- 
eralization replacing e ; with <j>i s , where s = 1, TV. [18] 
Taking the mean-field ansatz of ipi — p, iip r i = ip r , 
iQri = — Qr, an d i\ — A where Aj introduces the ro- 
tor constraint ^2 S \4*i S \ 2 = 1, we obtain the free energy 
functional 

Fmf = 



2 \ 

~~fj E E Mi^n - E k+ ) - — E ln( ^« ~ Ek -) 

ka Un ka W, 

+ \ E E ln (-7jK + vA 2 + v4 + a) 

ks V n 

+ ^2( Dx y - jjItVpt -p + y] +p + (fi r + q r 

k 

-A + /*[!-<*]). (10) 



Here the renormalized fermion spectrum is 
gi ven by E k ± = [ ^ zftzrt+toizgd ] ± 

^te^ P + V 2 , and 5 is hole con- 
centration for the conduction band. u> n (u n ) is the 
Matzubara frequency for fermions (bosons). 

Minimizing the free energy Eq. (10) with respect to 
A, x, y, q r , ip r , p, and /i, we obtain the self-consistent 
mean-field equations. Performing the Matzubara fre- 
quency summations and momentum integrals with ^2 k — 



■5 



Td J D d ^ e (constant density of states[5, 18]), we find 



(V/D, U/D) c is determined by 

V 2 r, f sinltan-^^ffn)]-! 



2U 
~~D 



2(L>/3) 2 



In 



sin [tan" 



''(tI 



■{ 



1 



l 



sin[tan- 1 ( £>/ |~ /f '- )] - 1 
sinftan-H^^)]- 



_ sin [ tan -i ( ^+^)] 



^U(\ + Dy)-^U(X + yeg) 
Dy 



co = 



;(#-*)■ 



.T = 



(2A - Dy) x /U(X + Dy) - (2A - yeg) x /U(X + yeg) 



-ln{^ 

^ — si 

-{ 



sin [ tan -l(^+^)] 
1 



sin [ tan -l(^:)] + l 



sin[tan- 1 (^±^)] + 



4V L> 2 > 

{= 1 



3(W 

V" 2 r. f sin[tan- 1 ( xe ^ r )] - 1 



'■) 2D 2 x 2 
^iN""^^)] 7 ! -sinltan-^^ftH 



sin [ tan -l(^)] 

1 



In 



r sin 
I- — si 



l} 



(5 = 



fJ-rV 

(D/3) 2 
V 



1 



cos[tan H ^V^ )] 



sin [ tan -i (±^)] + l 

1 ( Dx+H r 



D/3 

D 
3 



cos [tan 1 ( D ^y tJ ' r 
1 



cosftan-^^^)] 



3 V 



V 2 



[[tan" 1 ^ 



2V 



2V 



cos[tan-H^^)] 



(12) 



+ 



sin [t an -i(£S|^:)] + l _ si n[tan- 1 (^|±^)] + 
fj, r V 



L> 2 x 2 



cosltan-^^^)] costtan-^^i 



i / y 2 \ 

£V = 1" M + V ), 

X \ €J — Or- — Wr- / 



2 

— (</V - p + g r ) - 1 
V 



Actually, we find (0.112,0.044) at 5 = 0.800 and e d /D = 
-0.001, thus (V/U) c = 2.545. 

This quantum phase transition results from gapping 
of density fluctuations of localized fermions in the pres- 
ence of strong hybridization with renormalized conduc- 
tion fermions, thus differing from the Kondo breakdown 
transition[4] in the ALM. In V/U > (V/U) c collective 
density fluctuations become softened, and the sa ddle- 
point equation of A is replaced with 1 = yj2U/Dy + Z, 



where Z = \ (e l 



is the condensation amplitude. Thus, 



Dx 



cosftan-^^l^)] cosltan-^^t^)] 



|(^-p) = (l + |) = (l-g), 



fj, r = fi + p + e d - q r 



the coherence temperature is given by Tfl ~ (Dy — 
2U) 1 / 2 near the heavy-fermion Mott critical point, and 
below this temperature the valance-fluctuation-induced 
heavy-fermion phase arises. 

(H) V. HEAVY-FERMION SPIN LIQUID TO 

HEAVY-FERMION FERMI LIQUID QUANTUM 
TRANSITION: ITS CRITICAL FIELD THEORY 
AND NON-FERMI LIQUID PHYSICS 



Condensation of the density-fluctuation bosons occurs 
when their excitation gap closes or its effective chemical 
potential vanishes, i.e., ip r = 0. This causes p = 0, re- 
sulting in eg = —D (e^ = D). We obtain A = Dy at the 
quantum critical point, thus find the relation Dy/U = 2, 
completely consistent with the previous qualitative anal- 
ysis for the rotor model. [15] The hopping parameter for 
renormalized conduction fermions is x — 1/3 at the quan- 
tum critical point. As a result, the quantum critical point 



To investigate low energy physics near the heavy- 
fermion spin liquid to Fermi liquid quantum criti- 
cal point, we derive an effective field theory from 
Eq. (8). From the fermion mean- field Lagrangian 



E 



No- 



where 1p k<T 



Xka+(d T +E k+ )Xka+- 
= U k Xka+ + VkXka- 



-Xka-(d T +E k -)Xk<j- 

and t]ka = v k Xka+ 



UkXka- With life = V/ 



E k + - ef ) 2 + V 2 and v k 



-(E k+ - <%)/yJ (E k+ - ef) 2 + V 2 , we find its critical 
field theory C x = xL(d T - Hc)Xra + 2^H(V - 



ia r )xr 



Here the renormalized fermion field Xro- rep- 



resents Xicr- in the continuum limit, and the effec- 
tive band mass is given by m x — (d 2 Ek-/dk 2 ^J = 



with the "bare" band mass 
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m e ~ (tx)~ 1 , renormalized localized level e^r = £d + 
and critical chemical potential fi c . Note that m x diverges 
in the limit of V — > oo. 

The effective boson Lagrangian can be written with 
an electromagnetic field A as = [(d T + i<p r )(f>l s ][(d T — 
i(p r )(t>rs] + |(V - ia r - iA)4> rs \ 2 + TO 2 |</> rs | 2 + ^\4>rs\ 2 , 
where the rotor field e t6r is replaced with </> rs , and 
the rotor constraint is softened via introduction of lo- 
cal interactions u^. is an effective mass, given by 
to^ = (V/U) c — (V/U). An important point is that the 
renormalized boson chemical potential ip r becomes zero 
at the quantum critical point (to 2 = 0), thus the linear 
time-derivative term vanishes. 

We find the critical field theory at the heavy-fermion 
spin liquid quantum critical point 



/ 



drd 2 x 



1 



+ \d T <p rs \ 2 + |(V - iar - iA)</» rs | 2 + -^\<f> rs 
+Se ff [a r }, (13) 

where the critical gauge action is given by S e ff[a r ] — 
5 J2 q ,uj n ("fF 1 -^- + fqj (Sij - ^^CLirCLjr with the damp- 
ing strength jp = kp/n and boson flavor number 
TV. [18] The first term represents dissipative dynamics in 
gauge fluctuations due to particle-hole excitations of \ r(T 
fermions near the Fermi surface, and the second term 
arises from critical boson fluctuations at the quantum 
critical point. [18] The Maxwell gauge action is omitted 
since it is irrelevant in the renormalization group sense. 
Note that the time component of the gauge field mediates 
a local interaction, thus safely ignored in the low energy 
limit. [19] 

In the random phase approximation one finds the 
following expression for the free energy F/V = 



21 



tan 



ImD(ij,w) 
RcD(q,w) 



, where the gauge 



kernel is given by D(q,u>) = ijF^ + -f-q^ in real 

frequency at the quantum critical point. This free en- 
ergy leads to the divergent specific heat coefficient 7 = 
Cy/T oc — InT. The dc conductivity is given by the 
Ioffe-Larkin combination-rule a tot — a 4> G xl ( a 4> + a x)> as ~ 
suming gaussian gauge fluctuations. [20] In the one loop 
diagram of gauge-boson exchange one finds l/r t * ~ 
T 2 [18] and l/r? r ~ T[19], where 1/r* (l/rf r ) is the scat- 
tering rate of the Xra fermion (<j> rs boson) with tempera- 
ture T. As a result, the total dc conductivity is given by 
&tot ~ T in low temperatures, giving rise to the linear 
resistivity p to t ~ T at the quantum critical point. On the 
other hand, in V/U < (V/U) c but \V/U - (V/U) c \ 
corresponding to the heavy-fermion spin liquid state in 
the zero temperature limit, the gauge kernel shows the 

crossover behavior from D{q,u) = 

to D(q,u>) = (^—i"fF^ + with an effective inter- 

nal gauge-charge g as temperature goes down. 



Here 



the crossover energy scale would be an excitation gap 
of <j) rs , given by \<p r \. Accordingly, the specific heat coef- 
ficient exhibits the upturn behavior from 7 ~ — In T to 
7 ~ T -1 / 3 . But, the resistivity still behaves as p to t ~ T 
owing to p x ~ T 4 / 3 in the z = 3 critical field theory, [19] 
where z is the dynamical exponent. 



VI. DISCUSSION AND SUMMARY 

A. Stability of the heavy-fermion spin liquid 
against gauge fluctuations 

The stability of the heavy-fermion spin liquid state 
against gauge fluctuations is an important problem for 
this phase to realize in the finite-U ALM. It has been 
shown that the fermion-gauge action in Eq. (13), ob- 
tained via integrating out gapped boson excitations, has 
an infrared stable interacting fixed point with an effec- 
tive nonzero internal charge when the flavor number or 
density of fermions is sufficiently large. [21] Actually, the 
present author has studied that such a fixed point can 
arise when the fermion conductivity, given by its current- 
current correlation function, is sufficiently large. [22] No- 
tice that the fermion conductivity is associated with the 
density of fermions. At this fixed point the critical field 
theory is characterized by the dynamical exponent z = 3, 
as discussed before. 

The problem is the stability of this fixed point against 
instanton excitations resulting from compact gauge fluc- 
tuations. A similar situation appears in the compact 
QED3, where Dirac fermions interact via compact U(l) 
gauge fluctuations. The infrared interacting fixed point 
in the QED 3 has been shown to be stable against in- 
stanton fluctuations since the scaling dimension of an in- 
stanton insertion operator is proportional to the flavor 
number of massless Dirac fermions, thus irrelevant in the 
large flavor limit. [23] Following the similar strategy with 
this case, the present author has shown that the scaling 
dimension of the instanton operator is proportional to the 
fermion conductivity analogous to the flavor number of 
Dirac fermions, thus irrelevant in the large conductivity 
limit. [22] This seems to be natural because the fermion 
conductivity is associated with screening of gauge inter- 
actions. 

In the heavy-fermion spin liquid phase the fermion 
conductivity may not be sufficiently large because the 
strong hybridization makes the hybridized band flat. In- 
tuitively, the existence of the heavy-fermion spin liquid 
phase is not clear because the strong hybridization cou- 
pling is necessary for its existence, but such hybridiza- 
tion prohibits this phase from being stable against gauge 
fluctuations. Since its stability depends on the fermion 
conductivity, [22, 24] to determine its presence with gauge 
interactions is beyond the scope of the present study. 
However, its existence is allowed in principle. 

The stability of this phase against antiferromagnetic 
ordering is confirmed since the Fermi-surface nesting does 
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not exist in the heavy-fermion band. This not only justi- 
fies our neglect of spin fluctuations, but also propose new 
mechanism for the existence of the spin liquid phase. 



B. A variational wave function 

The above discussion motivates us to find the heavy- 
fermion spin liquid phase from the microscopic model 
itself, here the finite-U ALM. In this respect we propose 
a variational wave function for its existence. Since we 
start from the conduction-valance hybridized state, the 
following wave function is naturally proposed 



|*> = iM4, + «(*)<*L)|o>, 



(14) 



where |0) is the vacuum state and a(k) is the varia- 
tional parameter given by a(k) = V/{(e c k — \i — £d)/2 + 
y/(e% — Ji — £d) 2 /4 + V 2 ] in U = 0. Turning on on-site 
interactions to localized electrons, the variational ground 
state is proposed to be [25] 

I* > =Ppg\&) 

= 11,(1 - «44)rW( c L + a(A;)4 CT )|0 >(15) 

where Ppg — 11^(1 — nnf^nf^) is the partial Gutzwiller 
projector, suppressing charge-fluctuation effects, with an 
interaction-dependent k determined self-consistently as a 
function of C/. [26] 

Since this ground-state wave function captures both 
relevant interactions, V and U appropriately, physics in- 
cluded in Eq. (15) is expected to be essentially imposed 
in the slave-rotor effective gauge theory [Eq. (8)]. When 
gauge fluctuations are ignored as the saddle-point ap- 
proximation in Eq. (8), the |<j>) state is exactly recovered 
from the fermion Lagrangian Lp = L^, + L v + Ly because 
Lf gives rise to the same a(k) in |<E>). Note that the gauge 
transformation a a = e~ l6,i ?/w is an essential procedure 
for obtaining this strong hybridized dynamics. On the 
other hand, the partial Gutzwiller projection operator 
Ppg is simulated in the slave-rotor representation, Lg 
taking into account on-site interactions of localized elec- 
trons appropriately. This approach is parallel to that in 
the doped Mott insulator problem, where the Gutzwiller 
projected BCS wave function, i.e., the RVB state had 
been proposed, [27] and such a variational ground state 
was simulated in the context of the gauge theory. [1] It will 
be an interesting project to check whether this ground 
state wave function can result in both the fractionalized 
Fermi liquid and heavy-fermion spin liquid. 



C. Summary 

In this paper we have studied the finite-U ALM in the 
strong hybridization limit. The conventional treatment, 
integrating out conduction electrons to obtain an effec- 
tive Lagrangian in terms of localized fermions and col- 



U/D 




V/D 

FIG. 1: Schematic phase diagram of the finite-U Anderson lat- 
tice model without antiferromagnetism at zero temperature: 
HF-FL (heavy-fermion Fermi liquid), HF-SL (heavy-fermion 
spin liquid), and F-SL (fractionalized Fermi liquid) 



lective density-fluctuation bosons, describes the heavy- 
fermion quantum transition from the heavy-fermion 
Fermi liquid with a large Fermi surface to the fraction- 
alized Fermi liquid with a small Fermi surface, where 
the abrupt volume change of the Fermi surface is in- 
volved. On the other hand, the strong hybridization 
approach shows the heavy-fermion spin liquid to heavy- 
fermion Fermi liquid transition, where this type of Mott 
transition differs from the orbital-selective one via Kondo 
breakdown. 

We propose a schematic phase diagram of the finite-U 
ALM in Fig. 1, where HF-FL, HF-SL, and F-SL repre- 
sent the heavy-fermion Fermi liquid, heavy-fermion spin 
liquid, and fractionalized Fermi liquid, respectively. It is 
important to notice that this model has two indepen- 
dent parameters scaled by the conduction bandwidth, 
(V/D, U/D). In the small U/D limit the hybridization 
coupling is relevant to form a hybridized band, thus the 
heavy-fermion Fermi liquid arises. In the large U/D 
limit the Kondo breakdown transition has been shown 
to occur in the KLM or infinitc-U ALM, thus there is 
a critical V/D separating the fractionalized Fermi liq- 
uid ((c^riu) = 0) from the heavy-fermion Fermi liquid 
((cj-r/o-) 7^ 0). On the other hand, in the large V/D limit 
the heavy-fermion band is formed first ((^J.r/ CT ) ^ 0), 
and increasing U/D (marked by the arrow line) is ex- 
pected to cause the Mott transition ((e~ 10 ) = — > 
(c^rjcr) s=s (e~ te )(ip'lri cr ) = 0) in this heavy-fermion band. 

We proposed the variational ground-state wave func- 
tion for the finite-U ALM since the stability of the heavy- 
fermion spin liquid state against compact gauge fluctua- 
tions cannot be fully confirmed in the effective field the- 
ory approach, thus it is necessary to check its existence 
from the microscopic model. It will be interesting to ob- 
serve such a spin liquid state in the strong hybridization 
and large repulsion limits, where the heavy-fermion band 
still exists, but this corresponds to an insulator owing to 
charge gap. 

We close this paper discussing the possibility of 
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valance- fluctuation-induced superconductivity. [28] Re- 
cently, we have generalized the U(l) slave-rotor formu- 
lation of the Hubbard model into the SU(2) one, allow- 
ing not only local density fluctuations but also pairing 
excitations. [5] Both collective charge fluctuations form 
an SU(2) slave-rotor matrix field, where its off-diagonal 
components are associated with superconductivity. We 



believe that the SU(2) slave-rotor decomposition is also 
available to the finite-U ALM. This superconducting 
mechanism may explain physics in PuCoGa$, where su- 
perconductivity occurs with HF physics at the same 
time. [9, 29] This interesting possibility is under inves- 
tigation. 
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